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Abstract-  In  this  paper,  a  new  method  for  processing 
wideband  radar  data  is  presented.  To  perform  the  full  degree 
of  freedom  wideband  processing,  3-D  space-time  adaptive 
processing  (STAP)  needs  to  be  implemented,  which  involves 
intense  computational  burden.  One  approach  in  this  case  is  to 
do  subband  STAP  processing  and  combine  these  outputs.  In 
this  paper,  instead  of  traditional  subband  processing,  incoming 
wideband  data  signal  is  modulated  by  multiple  carriers, 
combined,  and  filtered  prior  to  processing  using  narrowband 
STAP.  This  method  offers  a  significant  decrease  in 
computation  burden  compared  to  the  subband  method. 

I.  Introduction 

One  way  to  understand  the  wideband  data  is  to  consider  it 
as  a  collection  of  narrowband  data  scenes.  Recall  that  the 
narrowband  data  vector  Xj^{t)  for  an  N  sensor,  M  pulse 

array  has  the  form  [1] 

=  co,)s_{9,j,  CO,)  (1) 

i  i 

where  .  represents  the  clutter  scatter  return  from  the 
{U  jy^  patch  in  the  field  of  view,  A{6^.,(o^  the  array 
factor  at  frequency  co^  .  Here  ^(6>. represents  the 
spatio-temporal  steering  vector  at  frequency  co^  given  by 

^AOup  =  c»A®A0,j,  <»,)  (2) 

where  the  spatial  steering  vector  is  a{0,  co^ )  is  given  by 


In  (2)  -  (5),  refers  to  the  operating  wavelength,  the 
platform  velocity  along  the  line  of  the  array  and  r 
represents  the  pulse  repetition  interval.  With  Xj^{t)  in  (1) 

representing  the  narrowband  data,  and  x(0  the  desired 
wideband  data,  we  have 

At)  =  Yj^At)  ■ 

k=l 

With  c(t)  representing  the  received  wideband  clutter  data 
and  /(t)  representing  the  received  wideband  target  data, 
the  total  receiver  signal  in  (6)  can  be  written  as 

A0  =  /(0  +  £(0-  (7) 

Assume  that  the  target  is  moving  with  a  relative  velocity  F 
at  an  arrival  angle  of  (both  parameters  are  unknown),  the 

MN  X 1  target  signal  / (t)  have  the  form  of 

A(i)  =  /(f-0'-i)7;-(/c-i)r,)  (8) 

where  the  spatial  delay  is  given  by 

^  c 


A9,  co,)  =  X  (3) 

and  the  temporal  steering  vector  b(co^  ,  given  by 


and  the  temporal  delay  is  given  by 


T  = 


2VTs\nO^ 


(10) 


b{co,^ ,  CO,)  =  [l,  ,  •  •  • ,  J .  (4) 

Here  d  represents  the  interelement  spacing  distance  and 
Doppler  co^  is  given  by 


\I2 


Hence  the  target  return  has  the  form 

where  represents  the  m*  pulse  return  given  by 
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L=[fit-iM-\)T,),  (12) 

In  the  wideband  case,  the  optimum  processor  is  a  whitening 
filter  H(z)  followed  by  a  matched  filter  [2].  This  whitening 
processing  is  shown  in  Fig.  1.  The  white  noise  w{t)  in  Fig. 
1  represents  the  whitened  interference. 


Whitening 

Filter 

H(z) 


r 

^y{t)  =  g{t)+w{t) 


White 

noise 


Fig.  1  Whitening  of  wideband  data. 


Using  (7)  -  (12),  in  the  frequency  domain,  the  whitening 
filter  output  is  given  by 


Y{(D)  =  F{(D)n{e^") 


1 

^-jcoT, 

-jco{N-\)T, 


=  a{  0,  c6) 


e^"^^a{e,co) 


a{6,  c6) 


-\-w(co) 


b(y,  co)®a{0,  co)  =  ^(0,  V,  co) 

F(co)U(e^^)^(0,V,  co)  ^  w(cd)  =  g  +  w 
(13) 


where  g  represents  the  target  output  and  w  represents  the 

interference  output.  In  this  case,  the  matched  filter  is  given 
by  g* .  Thus  the  optimum  wideband  processor  can  be 

shown  as  in  Fig.  2. 


Z  =  c*Y(co) 


Fig.  2  Optimum  wideband  processor. 


As  a  first  approximation,  if  we  use  a  narrowband  like 
whitening  filter  with  a  constant  term,  then 

h(z)  =  r;''"  ,  (14) 

where 

R  (0}>o  (15) 


represents  the  interference/clutter  covariance  matrix.  In  that 
case,  the  output  of  the  optimum  wideband  processor  is 
given  by 


Z=gYico) 

=  {l\0,  V,  ®)R;’'^)(r;’'^X(®)) 

=  V,  cd)R-J  X(cd)  =  W\cd)  X(cd) 

where  W*(cd)  represents  the  optimum  wideband  STAP 
processor  and  it  is  given  by 

W(co)  =R-Ji(0,V,cd).  (17) 

Notice  that  the  optimum  wideband  STAP  processor  in  (17) 
has  the  same  form  as  in  the  narrowband  case.  However,  it 
is  a  frequency  sensitive  processor  and  it  is  difficult  to 
implement  at  all  frequencies  simultaneously.  When  the 
whitening  filter  H(z)  involves  delay-lines,  the  structure  in 
(17)  becomes  more  complex. 

In  summary,  the  phase  delays  in  wideband  STAP 
processor  become  frequency  sensitive  filters.  The  optimum 
wideband  processor  must  be  compensated  at  all  frequencies 
simultaneously.  In  practice,  subband  schemes  can  be  used 
for  wideband  processing.  However,  this  scheme  is 
suboptimal  since  narrowband  processing  is  done  on  each 
subband. 


II.  Subband  Approach 

In  the  subband  method,  the  received  wideband  data  is 
expressed  as  sum  of  multiple  narrowband  signals  as  in  (6). 
A  bank  of  band-pass  filters  spanning  the  total  signal 
bandwidth  is  used  to  divide  the  signal  into  sub-bands.  If  the 
bandwidth  of  the  bank’s  filters  is  small  enough  then  the  sub¬ 
bands  approximate  narrowband  signals  and  can  be 
processed  using  narrowband  STAP. 

A  typical  filter  bank  is  shown  in  Fig.  3  which  uses  least 
square  approximate  FIR  low  pass  filters  modulated  to  the 
center  frequency  of  each  frequency  band.  By  modulating  a 
low  pass  filter  to  a  desired  frequency,  the  passband  of  the 
filter  includes  desired  frequencies. 

The  sampling  frequency  for  the  filter  shown  in  Fig.  3  is 
635MHz  with  251  taps  for  each  filter.  A  total  of  10  filters 
within  the  filter  bank  spans  80MHz  bandwidth  is  used.  Each 
filter  has  a  3dB  bandwidth  of  8MHz.  The  center  operating 
frequency  of  the  array  considered  here  is  at  435MHz  and 
80MHz  bandwidth  corresponds  to  18.4%  Fig.  4  shows  the 
structure  of  the  subband  method  that  use  K  beamformers. 
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Fig.  3  Twenty  sub-band  least  squares  FIR  filter  bank  ranging  from 
395MHz  to  475MHz. 
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Fig.  6  shows  the  average  of  the  10  subbands.  A  uniform 
array  with  14  sensors  and  16  pulses  is  considered  here.  The 
array  interelement  spacing  is  selected  to  be  half  wavelength 
at  center  frequency.  The  clutter  to  noise  ratio  and  target  to 
noise  ratio  are  set  at  40  dB  and  0  dB  respectively. 


Fig.  6  Average  Power  Spectrum  for  10  sub-bands  processed 
with  sample  matrix  inversion  with  diagonal  loading  and 
subarray/subpulse  smoothing  (SMIDLSASPFB). 
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(a)  Top  View 


(b)  Side  View 


III.  Modulation-Sum-Filter  Approach 

This  modulation-sum-filter  method  presented  here  for 
combining  wideband  data  is  in  a  sense  equivalent  to  the 
focusing  approach,  although  unlike  the  focusing  method, 
there  is  no  need  to  estimate  the  focusing  matrices  [4]. 
Instead,  the  multiple  modulations  are  used  to  focus  the 
wideband  data  into  a  narrow  band  region  where  they  are 
summed  and  processed  using  narrowband  STAP.  This 
procedure  is  illustrated  in  Fig.  7  below. 


Fig.  4  Subband  approach. 


Fig.  5  shows  the  4  frequency  band  outputs  of  the  subband 
approach  when  10  subbands  are  used.  The  sample  matrix 
inversion  with  diagonal  loading  and  subarray/subpulse 
smoothing  (SMIDLSASPFB)  technique  is  used  here  [3]. 


(c)  Subband  7,  Freq.  =  443  MHz 


(d)  Subband  10,  Freq.  =  467  MHz 


Fig.  5  Angle-Doppler  output  for  different  subbands. 


A'(m) 


Fig.  7  Illustration  of  subbands  being  modulated  to  a  common 
center  frequency 


By  aligning  the  signals  being  added  prior  to  band  pass 
filtering,  only  one  band  pass  filter  is  required.  The  signals 
are  aligned  by  modulating  each  subband  to  the  center 

frequency  as  shown  in  Fig.  8  where  Xi^co)  is  the  Fourier 

transform  of  x(^)  ,  is  the  Fourier  Transform  of 

modulated  by  ^  represents  the 

bandpass  filter. 

The  sum  of  the  modulated  signals  can  be  expressed  as  the 
signal  multiplied  by  a  sum  of  modulating  functions  shifting 
each  subband  to  the  common  center 

<F>^X(^CD- =  (1^) 
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This  combined  signal  is  finally  bandpass  filtered,  effectively 
resulting  in  a  single  subband  containing  the  averaged 
information  of  various  subbands.  Thus 

Z((o)  =  H,,{(o)Y(a)  (19) 

represents  the  filtered  output.  In  the  final  step,  traditional 
narrowband  STAP  processing  can  be  applied  to  the  filtered 
final  output  in  (19). 


Fig.  8  Illustration  of  modulated  versions  of  the  same  signal. 


In  practice,  the  total  number  of  carrier  frequencies  being 
used  and  the  bandwidth  of  the  bandpass  filter  are  free 
design  parameters.  In  what  follows,  a  set  of  simulation 
results  are  presented  to  exhibit  the  effects  of  these  two 
parameters. 


IV.  Simulation  Results 

Using  the  Modulation-Sum-Filter  approach,  a  single 
target  at  0^  in  the  filtered  data  generates  multiple  target 

vectors  corresponding  to  frequency  Or 

equivalently,  processing  the  data  at  co^  generates  multiple 
targets  at  0^,  0^^  where 


bandpass  filter  with  8MHz  bandwidth.  The  sample  matrix 
inversion  with  diagonal  loading  and  subarray/subpulse 
smoothing  technique  using  10  samples  is  shown  here. 
Observe  that  the  wideband  target  is  clearly  identified. 
However,  the  Angle-Doppler  spread  occurs  as  shown  here 
where  a  single  target  with  velocity  =  40m/s  appears  as  an 
“extended  targef  ’  with  different  velocities. 


Doppler  spread 


(a)  Top  View  (b)  Side  View 

Fig.  9  Angle-Doppler  spreading.  Data  is  modulated  by  10 
carriers  to  435MHz.  Combined  data  is  filtered  with  a 
bandpass  filter  with  8MHz  bandwidth. 


Fig.  10  shows  the  Angle-Doppler  output  using  the 
Modulation- Sum-Filter  approach  using  20  modulations. 
The  combined  data  is  filtered  using  a  bandpass  filter  with 
4MHz  bandwidth.  Once  again,  the  target  is  clearly 
identified  and  the  Doppler  spread  present  there  is  visible. 

On  comparing  Fig.  10  with  Fig.  9,  the  sidelobe  level 
using  20  modulations  is  lower  than  that  using  10 
modulations.  Thus  the  performance  of  the  processing 
output  using  20  modulations  is  superior  to  that  using  10 
modulations.  In  practice,  the  total  number  of  carrier 
frequencies  being  used  and  the  bandwidth  of  the  bandpass 
filter  need  to  be  selected  carefully  to  obtain  optimum 
performance. 


Doppler  spread 


(a)  Top  View  (b)  Side  View 

Fig.  10  Angle-Doppler  spreading.  Data  is  modulated  by  20 
carriers  to  435MHz.  Combined  data  is  filtered  with  a 
bandpass  filter  with  4MHz  bandwidth. 


co^  sin  0^=  sin  0j^,  ^  =  2,  3,  •  •  • ,  T .  (20) 

As  a  result,  “Angle-Doppler  spread”  occurs  in  STAP  output 
spectrum  processed  at  a  single  frequency. 

Fig.  9  shows  the  Angle-Doppler  output  using  the 
Modulation- Sum-Filter  approach.  In  Fig.  9,  ten  carriers  are 
used  for  modulating  the  data  to  the  center  frequency  435 
MHz.  The  combined  data  is  then  filtered  with  a  single 


V.  Conclusions 

This  paper  presents  a  new  method  for  processing 
wideband  radar  data.  Instead  of  traditional  subband 
processing,  incoming  wideband  data  signal  is  modulated  by 
multiple  carriers,  combined,  and  filtered  prior  to  processing 
using  narrowband  STAP.  This  method  offers  a  significant 
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decrease  in  computations  compared  to  the  subband  method. 
The  detection  performance  is  affected  by  free  parameters 
such  as  the  number  of  modulations  used  and  the  bandwidth 
of  the  filter.  Using  the  Modulation- Sum-Filter  approach,  a 
single  target  in  the  filtered  data  generates  multiple  targets 
corresponding  to  different  frequencies.  As  a  result,  “Angle- 
Doppler  spread”  occurs  in  STAP  output  spectrum  processed 
at  a  single  frequency.  Methods  to  align  the  angle-Doppler 
spectrum  need  further  study. 
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